Characterization of Disorder in Semiconductors via Single— Photon Interferometry 
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The method of angular photonic correlations of spontaneous emission is introduced as an ex- 
perimental, purely optical scheme to characterize disorder in semiconductor nanostructures. The 
theoretical expression for the angular correlations is derived and numerically evaluated for a model 
system. The results demonstrate how the proposed experimental method yields direct information 
about the spatial distribution of the relevant states and thus on the disorder present in the system. 

PACS numbers: 78.55.-m, 42.50.-p, 71.35.-y, 78.30.Ly 



Disorder in solids is a topic of increasing interest in 
solid-state research in the last decades. Intriguing and 
widely celebrated phenomena likeQuantum-Hall-Effect, 
photon echoes, and Anderson localization rely on the ex- 
istence of disorder. However, the significance of disorder 
is not purely academic. Modern devices are fabricated 
on smaller and smaller length scales. Systems with re- 
duced dimensionality can be realized by semiconductor 
heterostructures. Concomitantly, disorder due to inter- 
face roughness and/or compositional fluctuations plays 
a crucial role in the performance of these nanostruc- 
tures. Clearly, for a detailed understanding of the in- 
fluence of disorder on the investigated phenomena, it is 
important to know to what degree the electronic states 
deviate from the Bloch-like character which governs the 
electronic properties in perfectly periodic structures. In 
particular, localized electronic states can occur which 
strongly alter both the optical and transport properties 
of a given system [H, HI 0, 0, IE IS 0- ^ is, therefore, 
important to have a reliable method that allows the ex- 
perimentalist to quantitatively characterize the disorder 
in a given structure without its destruction. Unfortu- 
nately, in this respect the present situation is not very 
satisfactory. 

For example, the lattice disorder in alloy semiconduc- 
tors can be determined in princi ple by spatially resolved 
atomic-level microscopies [1, lil. lldL UlL Il2j |. but its in- 
fluence on ejr. optical properties is still not well under- 
stood [1, 0, 0] . Alternatively, one can obtain information 
about the connection of disorder and optical properties 
directly within the resolution limit set by the wavelength 
of light. Examples include micro- (nano)-lumincscence 
[HI EH i near-field optical ima g i ng fla. jltij . and resonant 
Rayleigh scattering (RRS) [l7|, [Tj| In order to de- 
duce the spatial dependence of the disorder potential, the 
first two methods require a spatial scanning of the sam- 
ple. When the angular dependence of the RRS intensity 
is measured, one can determine the energy statistics of 
emitters (HI, as well as the localization length scales of 
the optically active excitonic states [l7, HI, ljj]. However, 



RRS can deduce the exact disorder potential and posi- 
tion of emitters only indirectly via an inverse-scattering 
problem. 

In this Letter, we suggest a quantum-optical interfer- 
ometric method which provides the disorder potential 
directly. Our scheme is based on the photonic correla- 
tions that are determined from the stationary sponta- 
neous emission by collecting two different emission direc- 
tions into a common detector. Recent experiments have 
shown that strong single-photon interferences can be ob- 
served in this setup [20( ■ We demonstrate that the qual- 
itative features of these interferences can be utilized to 
reconstruct the disorder landscape of the light emitters 
for length scales varying between the used wavelength 
and the illumination spot size. We also show that the re- 
construction of disorder potential can be performed solely 
on the basis of experimental data and that the method 
is robust to experimental inaccuracies. In contrast with 
micro- (nano)-luminescence fl3L 14 1 and near- field opti- 



cal imaging [lEI, \±j^, our scheme does not need spatial 
scanning. 

To present the idea behind our proposed quantum- 
optical method, we need to quantize the light field and 
compute the correlations among the emitted photons in- 
cluding the disorder as well as the relevant many-body 
interactions among the electron-hole excitations in the 
semiconductor heterostructure. For this purpose, and in 
order to keep our analysis as transparent as possible, we 
choose a simple model based on a one-dimensional tight- 
binding description. For a chain of N sites with periodic 
boundary conditions, every site, if isolated from the oth- 
ers, supports one electron and one hole state with the 
energies e| and , respectively. The index j refers to the 

position of the site and the energies e^' vary from site- 
to-site due to disorder. The corresponding carrier system 
is described by the material Hamiltonian 
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e,- and h[,hj denote the electron and hole cre- 



ation and annihilation operators at a site j . The constant 
J e ' h defines the bandstructure via the nearest neighbor 
coupling indicated by the symbol < jl >. We concentrate 
here on low-density conditions such that we only describe 
the attractive Coulomb interaction of electron-hole pairs, 
given by the last term in Eq. ([1]) . The model parameters 
are chosen to reproduce the effective mass and the ex- 
citon binding energy of typical GaAs quantum-well sys- 
tems In particular, iV sites are positioned evenly at 
distances a with J e = —5.737 meV and J h = 0.835 meV. 
The regularized Coulomb matrix element is given by 
Vji = Uo/(\j — l\a + ao) with Uq = 7 meV and ao = 0.5 a. 
This yields an exciton binding energy of Ef, = 7.5 meV. 

The quantized light field can be expressed via 
a mode expansion in terms of plane-wave modes 

momenta q 21, \22 
The light-matter 



photon 
tion is 



e iqr with 
operators /3 q . 
given by H h 



and 
interac- 

E q ^(B q S q + l/2) - 
3 (irjUqjBqejhj + h.c), where the absorption of a 
photon creates an electron-hole pair at the site j since 
only direct transitions are allowed. This term is propor- 
tional to r) = hq£ with the dipole-matrix element /j,q 
and the vacuum-field amplitude £. The light-matter in- 
teraction also depends on the mode function u qj at the 
position of the site, Rj . The sites can be assummed to 
be equally spaced since their size is much smaller than 
the optical wavelength. Consequently, we may choose 
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such that all disorder effects enter 



1,3 



via the energies e -' . The corresponding u, 
contains the in-plane photon momentum in the direction 
e,|. 

Under incoherent conditions, all coherent quantities, 
(hjej), (eljhj) and (-B q ) , (B q ), vanish. In this situa- 
tion, the system emits light only spontaneously leading 
to non-vanishing (B^ Bq>). For steady-state conditions, 
the photon-flux d/dt{B^Bq) defines the photolumines- 
cence (PL) intensity while d/dt{B^Bqi) determines the 
quantum-optical correlations between two different emis- 
sion directions. By using a suitable optical arrangement, 
one can collect light in the directions q and q' ^ q to 
a common detector. With this interferometric setup, 
E/q,q' = d/dt{B\ t-Bq') defines the contrast of the single 
photon interferences (20| for steady-state emission con- 
ditions. This quantity follows from 
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showing that J7 q , q ' is coupled to the photon-assisted po- 
larization n q ; = (B^hjei) and its dynamics 
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where - ef + t) - £ m (W£ + *W£)- The sec- 
ond line in Eq. (3) contains contributions due to Coulomb 
scattering and stimulated emission. The stimulated term 
can be omitted for systems without a cavity [22^ . As 
another simplification, we introduce a constant dephas- 
ing ftjjn q ; | scatt = — 7il q ; . When the carrier system is 
in a quasi-equilibrium, Sji = Y, m u i,™Ju Ci j(e\ n hl n h : jei) 
acts as a constant spontaneous-emission source and the 
Coulomb term, VII, leads to excitonic resonances in the 
emission. At the lowest level, Sji is given by its Hartree- 
Fock approximation Sji = djifjfj 1 containing the excited 
state occupations, /? = (ejeA and /j 1 = (h/-hj). 

In order to understand how disorder effects enter in 
^q,q': we nrs t solve Eqs. ([I])-© for a simplified case by 



setting all carrier-interaction terms to zero J 
0. This yields a steady-state expression 
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where ej = e| + e'j. We have also assumed that |q| = |q'| 
differ only by the emission direction giving us = c |q| and 
Aq = <7|| — q',, . The normalization j\f is chosen such 
that the correct 100% interference contrast is observed 
for q = q'. This form clearly shows that even the spon- 
taneous emission contains a phase, e lAqRj , influenced by 
the position of the emitters. This phase can always be 
observed as single-photon interferences [13] , which is the 
fundamental principle behind the proposed interferome- 
try. 

In the continuum limit a — > 0, Eq. ([4]) reduces to 



U hw (Aq) =jV J 
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lAqx dx . 
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where S(x) defines the level of excitation. We now notice 
that Fourier transformation of Eq. ([5|) with respect to Aq 
produces a peak at e(x) = has; i.e. the position of the 
emitters can be determined. Thus wc define 
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where Ago = 2lu/c = Att/Xq is the largest in-plane mo- 
mentum which is accessible in principle. This allows us 
to reconstruct the disorder potential up to the limits set 
by Ago via 



U(x) 



E2 

has I Uhui (x) I dhus/ / I Uhu (x) \ dims . (7) 
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The integration limits are set by the relevant spectral 
range E\ < hes < E 2 - Equations ©-Q suggest a 
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FIG. 1: a) Correlation function Unui=E Q (Aq) as a function of 
Aq/ Aqo for noninteracting TLS. The inset shows PL for the 
ordered (line) and the disordered (shaded) case, b) Contour 
plot of the spatially dependent correlation function Uhu>(x). 
Solid line is the input potential with C = 1.2Ao. c) The same 
for the fully interacting system, C = 0.2Ao. Eb is the sum of 
Eq and the binding energy, d) Dashed line: input potential 
with L = 0.25Ao, based on the same initial distribution of site 
energies for noninteracting TLS as b). Solid line: U(x), c.p. 
with solid line in b). 

general interferometric scheme to determine the disor- 
der landscape. First, the interference contrast U%uj (Aq) 
has to be measured as function of energy for different 
emission directions. Second, a simple Fourier transform 
of Unu) (Aq) produces directly the disorder landscape of 
emitters. 

In our full numerical computations, we generate the 
energies e| and from a random box-like distribution 
of width W e and W h , respectively. We take W e /W h = 
\J e I 'J |, which guarantees that the disorder amplitudes 
scale with the effective masses. The ordered case is re- 
covered by setting W e = W h = 0, i.e., ef l = e e ' h . The 
length scale of the disorder potential C is defined by a 
Gaussian low-pass filter applied in Fourier space with a 
maximum spatial frequency cx C^ 1 . We use Ao = 800 nm, 
band gap energy E = 1.5 eV, 7 = 0.1 mcV, N = 2048, 
and a = 5 nm. 

To illustrate the single-photon interferometric scheme, 
we first consider the non-interacting case and set the ex- 
citation S(x) to a constant value within a 10 /im exci- 
tation spot. The inset in Fig. [IJi shows the photolumi- 
nescence for W e = 0.5 meV (shaded area) and the or- 
dered case (solid line); clear inhomogencous broadening 
is observed. The interference contrast, Utuj (Ac?) is pre- 
sented in Fig. QJi for fiu = Eq as function of Aq. The 
central peak at Aq = produces a 100% interference 
constrast while U drops quickly for elevated \Aq\. The 
peak structure at small |Aq| values contains direct infor- 
mation about the disorder. The width of both central 
and side peaks is 2ir/S where S is the size of the spot. 



After Unu) (Aq) is defined for all relevant frequencies, we 
may construct Uhu>(x) via Eq. ((6|). We use the maximum 
optically allowed Aq (denoted by shaded area in Fig. Q}i) 
as the integration limits. The generated U%u (x) > i 
contour is shown in Fig. [Tb as shaded area, denoting the 
full width of half maximum of Uhu> {%), while the original 
potential is given by the solid line. Clearly the center 
of Uhu> (x) reproduces excellently the x dependence of the 
actual disorder landscape. For a fixed x, Uhuj(x) gives the 
homogeneous line for the /iPL spectrum of the emitter at 
position x. In the present paper, this width is deter- 
mined by the constant 7. In reality, microscopic carrier 
and phonon scattering can lead to a position dependent 
homogeneous line width which should be observable in 
the same way. The total inhomogeneous spectrum fol- 
lows as an integral over x. In Fig. [TJi we show results 
for a disorder potential with C = 0.25Ao- The long-range 
features of this potential are very well reproduced by our 
method. 

From an experimental point of view it is critical to 
know i) the range of q-vectors which can be detected 
and ii) how accurately the side peaks of U must be mea- 
sured with respect to the central peak. Concerning i), 
the experimental setup from Ref. [17| would allow for a 
range of Aq = N.A./n GaAs ^ 0.12Aq only (N.A. is the 
numerical aperture). Augmenting the setup by a solid 
state immersion lens with high index of refraction (for 
a discussion see [24|), and using two optimal positioned 
microscope objectives Aq ~ ±0.85Aqo can be achieved. 
This range is indicated by the shaded area in Fig. QJi. 
Concerning ii), we have solved Eq. (|5|) analytically for 
a sinusoidal e(x). In particular, we find that the sec- 
ond peak of Ufru (Aq) is located at Aq = 4n/C and the 
ratio r between the central peak and the second one is 
r sin = (v 7 ! + (7/AWO 2 + 7/AVF)- 2 , where AW is the 
standard deviation of the potential. Taking two sinu- 
soidal functions with equal amplitude but different spa- 
tial frequency yields values for r situated clearly below 
r s i n . For various random potentials, we find r is around 
r s ; n /2. This analysis shows that the ratio 7/ AW defines 
the accuracy requirements in an experiment. For exam- 
ple, if one wants to resolve disorder up to AW > 7/2, 
one has to measure the interference contrast with 10% 
accuracy. 

To verify the general validity of the scheme defined 
by Eqs. ((6j)— ((TJ) , we next investigate the fully interact- 
ing case by including the band-structure, disorder, and 
Coulomb interaction terms in Eq. ([3]) . We solve the eigen- 
states and eigenvalues eA of the homogeneous part 
of Eq. |3]) for steady-state conditions. This procedure 
yields [25[ 

x d i K E \ - M + 7 2 
for constant S. We evaluate Eq. ([5]) numerically for a 
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FIG. 2: Sensitivity of the reconstructed potential to various 
external parameters: a) sensitivity of measurement of the in- 
terference contrast, b) random noise, and c) range of momen- 
tum scan. Details are in the text. 

disorder realization with W = 5j. Since our main pur- 
pose is to show that the interaction contributions do not 
affect the reconstruction procedure I©-©, we can re- 
strict the calculations to a rather small system (N = 140, 
a = 10 nm) to obtain numerically feasible computations. 
Physically, this implies a small S = 1.4 ^m spotsize. Nu- 
merically, this leads to a relatively large discretization in 
Aq which forces us to use the Aq data taken from the en- 
tire Brillouin zone in the reconstruction procedure. Fig- 
ure [TJ; shows the obtained full U%u (x) > | in a contour 
plot (shaded area) and compares it with the actual poten- 
tial (solid line) [26(. Again, the center of Un u (x) matches 
well with the original disorder potential, which demon- 
strates the general applicability and the robustness of the 
proposed scheme. 

To determine the sensitivity of the reconstruction pro- 
cedure on several experimentally relevant inaccuracies, 
we apply Eq. j5]) to reconstruct the potential via Eq. (J7|) 
for different situations of the non-interacting case. In 
all frames of Fig. [5] the shaded area indicates the orig- 
inal potential with W = 57 and the reconstructed po- 
tential U(x) is shown as a solid line. In Fig. Uhuj(x) 
is constructed by omitting values of Uhu> (Aq) with in- 
terference contrast below 10%. In Fig. [2b, Ufr u (x) is 
constructed from Uk u (Aq) where we simply added a 5% 
random noise. When the momentum scan is limited to 
\Aq\ < Aqo/2 and |Ag| < 2Aqg (dashed line) we obtain 
the Uhui(x) results presented in Fig. [2J:. This is also ac- 
companied by a limitation of the spatial resolution of the 
resulting reconstructed potential. For a minimum Aq m i n 
given by experimental limitations, the spot size should 
be chosen accordingly, such that S = 27r/Aq m - m in order 
to avoid aliasing problems. 

All these cases indicate that the proposed photon- 
correlation interferometry produces the disorder poten- 
tial even when a reasonable amount of experimental in- 
accuracies are present. 

In conclusion, we suggest that the experimental 



method of single-photon interferometry can directly 
characterize the optically relevant disorder landscape in 
semiconductor heterostructures. The method is only lim- 
ited by the wavelength of light and the optical spotsize. 
We have also demonstrated a remarkable quality of the 
reconstructed potential even for weak disorder and added 
random noise. 
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